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The Coulomb drag in a system of two parallel layers is the result of electron- 
electron interactions between the layers. We have observed reproducible fluc- 
tuations of the drag, both as a function of magnetic field and electron concen- 
tration, which are a manifestation of quantum interference of electrons in the 
layers. At low temperatures the fluctuations exceed the average drag, giving 
rise to random changes of the sign of the drag. The fluctuations are found to be 
much larger than previously expected, and we propose a model which explains 
their enhancement by considering fluctuations of local electron properties. 

In conventional measurements of the resistance of a two-dimensional (2D) layer an electrical 
current is driven through the layer and the voltage drop along the layer is measured. In contrast, 
Coulomb drag studies are performed on two closely spaced but electrically isolated layers, 
where a current I± is driven through one of the layers (active layer) and the voltage drop V 2 is 
measured along the other (passive) layer (Fig. 1). The origin of this voltage is electron-electron 
(e-e) interaction between the layers, which creates a 'frictional' force that drags electrons in 
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the second layer. The ratio of this voltage to the driving current Rr> = —V2/I1 (the drag 
resistance) is a measure of e-e interaction between the layers. The measurement of Coulomb 
drag in systems of parallel layers was first proposed in Ref. (|7]|2]) and later realised in a number 
of experiments (TJ] 0131^0) (for a review see Ref. ©). As Coulomb drag originates from e-e 
interactions, it has become a sensitive tool for their study in many problems of contemporary 
condensed matter physics. For example, Coulomb drag has been used in the search for Bose- 
condensation of interlay er excitons ©, the metal-insulator transition in two-dimensional (2D) 
layers (|70l) . and Wigner crystal formation in quantum wires (|77|) . 

Electron-electron scattering, and the resulting momentum transfer between the layers, usu- 
ally creates a so-called 'positive' Coulomb drag, where electrons moving in the active layer drag 
electrons in the passive layer in the same direction. There are also some cases where unusual, 
'negative' Coulomb drag is observed: e.g. between 2D layers in the presence of a strong, quan- 
tising magnetic field (00); and between two dilute, one-dimensional wires where electrons are 
arranged into a Wigner crystal (|77l) . All previous studies of the Coulomb drag, however, refer 
to the macroscopic (average) drag resistance. Recently there have been theoretical predictions 
of the possibility to observe random fluctuations of the Coulomb drag (|72]|73l) . where the sign 
of the frictional force will change randomly from positive to negative when either the carrier 
concentration, n, or applied (very small) magnetic field, B, are varied. 

Drag fluctuations originate from the wave nature of electrons and the presence of disorder 
(impurities) in the layers. Electrons travel around each layer and interfere with each other, after 
collisions with impurities, over the characteristic area ~ I? , where L v is the coherence length 
(Fig. 1). This interference is very important for conductive properties of electron waves. For 
example, the interference pattern is changed when the phase of electron waves is varied by a 
small magnetic field, producing universal conductance fluctuations (UCF) seen in small samples 
with size L ~ L v . There is, however, a significant difference between UCF and the fluctuations 
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of the drag resistance. The former are only a small correction to the average value of the 
conductance: in our experiment the single-layer resistance fluctuates by ~ 200 mOhm around 
an average resistance of approximately 500 Ohm. In contrast, the drag fluctuations, although 
small in absolute magnitude (~ 20 mOhm) are able to change randomly, but reproducibly the 
sign of the Coulomb drag between positive and negative. Surprisingly, we have found that 
these fluctuations of the Coulomb drag, observed at temperatures below 1 K, are four orders of 
magnitude larger than predicted in Ref. ((72j) . 

Our explanation of the giant drag fluctuations takes into account that, unlike the UCF, the 
drag fluctuations are not only an interference but also fundamentally an interaction effect. In 
conventional drag structures the electron mean free path / is much larger than the separation d 
between the layers, and therefore large momentum transfers hq between electrons in the layers 
become essential. According to the quantum mechanical uncertainty principle, ArAg ~ 1, 
electrons interact over small distances Ar <C / when exchanging large values of momentum 
(Fig. 1). As a result the local properties of the layers, such as the local density of electron states 
(LDoS), become important in the interlayer e-e interaction. These local properties at the scale 
Ar <C I exhibit strong fluctuations (|7?|) that directly manifest themselves in the fluctuations of 
the Coulomb drag. 

The samples used in this work are AlGaAs-GaAs double-layer structures, in which the car- 
rier concentration of each layer can be independently controlled by gate voltage. The two GaAs 
quantum wells of the structure, 200 A in thickness, are separated by an Alo.33Gao.67 As layer 
of thickness 300 A. Each layer has a Hall-bar geometry, 60 pm in width and with a distance 
between the voltage probes of 60 /im (|75l) . 

Figure 2 shows the appearance of the fluctuations in the drag resistivity, pn, at low temper- 
atures. At higher temperatures, the drag resistance changes monotonically with both T and n: 
the insets to Fig. 2 show that pn increases with increasing temperature as T 2 and decreases 
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with increasing passive-layer carrier concentration as n b 2 , where b ~ —1.5. These results are 
consistent with existing experimental work on the average Coulomb drag (l?]|7oT). 

Figure 3A shows a zoomed-in view of the reproducible fluctuations as a function of n 2 . 
These fluctuations result in an alternating sign of the drag, which is demonstrated in the inset 
to Fig. 3 where the temperature dependence of the drag is shown at two different values of n 2 . 
The drag is seen first to decrease as the temperature is decreased, but then become either in- 
creasingly positive or increasingly negative, dependent upon n 2 . The reproducible fluctuations 
of the drag resistivity have also been observed as a function of magnetic field (Fig. 3B). For a 
fixed temperature, the magnitude of the drag fluctuations as a function of n 2 is roughly the same 
as that as a function of B. 

The theory of Ref. (|72)) calculates the variance of drag fluctuations in the so-called diffusive 
regime, I < d. In this case the drag is determined by global properties of the layers, aver- 
aged over a region Ar ^> /. The expected variance of drag fluctuations (at low T when the 
fluctuations exceed the average) in the diffusive regime is 



e 4 Et(L)t,p In nd 



where od ~ Pd/{piP2), and p\ and p 2 are the active and passive layer resistivities, respec- 
tively; Et(L) is the Thouless energy, Et(L) = hD/L 2 , D is the diffusion coefficient; t v is 
the decoherence time; k is the inverse screening length; A = 4.9 x 1CT 3 and g = h/(e 2 p) is 
the dimensionless conductivity of the layers. Using the parameters of our system, this expres- 
sion gives a variance of ~ 6 x 1CT 11 fiS 2 , which is approximately eight orders of magnitude 
smaller than the variance of the observed drag fluctuations. The fluctuations in pr, have been 
measured in two different samples, and their variance is seen to be similar in magnitude and 
T-dependence, confirming the discrepancy with the theoretical prediction (fi2|) . 

The expected fluctuations of the drag conductivity share the same origin as the UCF in 



4 



the conventional conductivity: coherent electron transport over in the layers prior to e-e 
interaction between the layers (Fig. 1). For this reason we have compared the drag fluctuations 
with the fluctuations seen in the single-layer resistivity of the same structure (Fig. 3B, inset), 
which have shown the usual behaviour (|771) . We estimate the expected variance of the single- 
layer conductance fluctuations using the relation (Ao^J) = (e 2 / h) 2 [Lt / L) 2 ', where Lt = 
J%D /ksT is the thermal length (17). This expression produces a value of 0.8 fiS 2 , which is in 
good agreement with the measured value of 0.6 /iS 2 . The typical 'period' of the drag fluctuations 
(the correlation field, AB C ) is similar to that of the UCF £75]), indicating that both depend upon 
the same L 9 and have the same quantum origin. 

To address the question of the discrepancy between the magnitude of drag fluctuations in 
theory (Ti2l) and our observations, we stress that the theoretical prediction for the variance, Eq.[H 
was obtained under the assumption of diffusive motion of interacting electrons, with small 
interlayer momentum transfers, q <C 1//. As the layers are separated by a distance d, the 
e-e interactions are screened at distances Ar > d. Therefore, in all regimes the maximum 
momentum transfers are limited by q < 1/d. In the diffusive regime, I < d, this relation also 
means that q < l/l, that is, interlayer e-e interactions occur at distances Ar > / and involve 
scattering by many impurities in the individual layers. In the opposite situation, I ^> d, the 
transferred momenta will include both small and large q- values: q < l/l and l/l < q < 1/d. 
We have seen that small q cannot explain the large fluctuations of the drag (12 ), and so argue that 
it is large momentum transfers with q > l/l which give rise to the observed effect. In this case 
the two electrons interact at a distance Ar that is smaller than the average impurity separation 
and, therefore, it is the local electron properties of the layers which determine e-e interaction. 
In Ref. CZU) it is shown that the fluctuations of the local properties are larger compared to those 
of the global properties that are responsible for the drag in the diffusive case. 

A theoretical expression for the drag conductivity is obtained by means of a Kubo formula 
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analysis (]7^1 [7P1 1201 1271) (detailed description in supporting text). For a qualitative estimate, 
three factors have to be taken into account: (i) the inter-layer matrix elements of the Coulomb 
interaction T>^\ (ii) the phase space (the number of electron states available for scattering); and 
(iii) the electron-hole (e-h) asymmetry in both layers. Point (iii) takes into account that in a 
quantum system the current is carried by both electron-like (above the Fermi surface) and hole- 
like (below the Fermi surface) excitations. If they were completely symmetric with respect to 
each other, then the current-carrying state of the active layer would have zero total momentum 
and thus no drag effect would be possible. The physical quantity that measures the degree of e-h 
asymmetry is the non-linear susceptibility T of the 2D layer. Theoretically, the drag conductivity 
is represented in terms of the non-linear susceptibilities of each layer and dynamically screened 
interlayer Coulomb interaction V i3 -.(u) as oc I dujVi2(uj)T 2 (uj)V 2 i(uj)Ti(uj) (indices 1 and 
2 correspond to the two layers) (17^11721) . The e-h asymmetry appears in Y as a derivative of 
the density of states v and the diffusion coefficient D: Y oc d (vD) /dfi, and it is this quantity 
that is responsible for the fact that drag fluctuations can exceed the average. As Dv ~ g and 
the typical energy of electrons is the Fermi energy, Ep, we have d(uD)/dfi ~ g/Ep for the 
average drag. The typical energy scale for the interfering electrons, however, is E T (L^,) (1771) . 
which is much smaller than E F and therefore a mesoscopic system has larger e-h asymmetry. 

Under the condition of large momentum transfer between the layers, fluctuations in Y are 
similar to the fluctuations of the LDoS, which can be estimated as 5u 2 ~ (z/ 2 ; / 1 g) In (max(L ¥ ,, L T )/l) 
(|7?1) . Also, the interaction in the ballistic regime can be assumed to be constant, Z>y « — 1/ vnd, 
as q is limited by q < 1/d. Finally, to average fluctuations of the drag over the sample with 
size L we should divide it into coherent patches of size that fluctuate independently and thus 
decrease the total variance: (Acr%) = (AaKL^iL^/L) 2 . If k B T > E T (L V ), fluctuations 
are further averaged on the scale of ~ k B T, and therefore the variance is suppressed by an 
additional factor of Et^L^) /ksT . Combining the above arguments we find 
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where iV is a numerical coefficient. 

Compared to the diffusive situation (Eq. Q]) the fluctuations described by our model are 
greatly enhanced. The difference between Eqs. [2] and Q] comes from the fact that in the ballistic 
regime electrons are not scattered by impurities between events of e-e scattering. Large momen- 
tum transfers correspond to short distances, and thus in the ballistic regime drag measurements 
explore the local (as opposed to averaged over the whole sample) non-linear susceptibility. This 
leads to the appearance of three extra factors in Eq. [2] (i) the factor I 4 / d A (which is also present 
in the average drag in the ballistic regime - see Ref. diffl) ): (ii) the phase space factor T/Et 
(which appears due to the fact that interaction parameters are now energy-independent); 
and (iii) the extra factor g 2 due to fluctuations of the local non-linear susceptibility. Local fluc- 
tuations are enhanced since the random quantity T is now averaged over a small part of the 
ensemble, allowing one to detect rare impurity configurations. 

Our model not only explains the large magnitude of the fluctuations, but also predicts a 
non-trivial temperature dependence of their magnitude. The latter comes from the change in 
the temperature dependence of L v (|22l) : at low temperatures, k B Tr/h <C 1, the usual result is 
Ltp oc T -1 / 2 , while for ksTr/h > 1 the temperature dependence changes to oc T _1 (|23l) . 
Consequently, the temperature dependence of the variance of the drag fluctuations is expected 
to change from T _1 at low T, to T -4 at high T. This temperature dependence is very different 
from the T-dependence of drag fluctuations in the diffusive regime, (Aerfj) oc T _1 . To test the 
prediction of Eq. [2l the T-dependence of (Acr^) has been analysed (Fig. 4). The variance is 
calculated in the limits of both the diffusive r v (solid line, r^ 1 oc T) and ballistic r v (dashed 
line, r^ 1 oc T 2 ), using iV ~ 1CT 4 . In fitting the drag variance we have found r v to agree 
with theory to within a factor of two (|75l) . which is typical of the agreement found in other 



experiments on determining r v (|2?1) . (The single-layer values of found from our analysis of 
the UCF agree with theory to within a factor of 1.5.) Thus, the temperature dependence of the 
observed drag fluctuations strongly supports the validity of our explanation. 

We have observed reproducible fluctuations of the Coulomb drag and demonstrated that they 
are an informative tool for studying wave properties of electrons in disordered materials, and 
the local properties in particular. Contrary to UCF which originate from quantum interference, 
fluctuations of drag result from an interplay of the interference and e-e interactions. More the- 
oretical and experimental work is required to study their manifestation in different situations. 
For instance, similarly to the previous extensive studies of the evolution of UCF with increasing 
magnetic field, such experiments can be performed on the fluctuations of drag. One of the im- 
portant developments in the field of Coulomb drag fluctuations can be their study in quantising 
magnetic fields, including the regimes of integer and fractional quantum Hall effects. 
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Fig. 1. Schematic showing the origin of the drag signal V 2 induced by the current li. 
The fluctuations of the drag arise from the interference of electron waves in each layer, 
before the two electrons take part in the interlayer interaction. 
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Fig. 2. Drag resistivity as a function of passive-layer carrier concentration for different 
temperatures: T = 5,4,3,2, 1,0.4, and 0.24 K, from top to bottom. Inset (A): p D as a 
function of T 2 . Inset (B): p D as a function of n 2 , with ni = 1.1 x 10 11 cm 2 ; dashed line 
is a n 2 1 5 fit. 
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Fig. 3. (A) Drag resistance measured at low temperatures as a function of passive- 
layer concentration; T = 1,0.4, and 0.24 K, from top to bottom. Inset: p D as a function 
of T for two values of n 2 denoted by the dotted lines in Fig. 3A; solid line is the expected 
T 2 dependence of the average drag. (B): p D as a function of B; T = 0.4, 0.35, and 0.24 
K, from top to bottom. (Graphs for higher T are vertically offset for clarity.) Single-layer 
concentration for each layer is 5.8 x 10 10 cm 2 . Inset: The UCF of the single-layer, with 
an average background resistance of 500 Ohm subtracted. 
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Fig. 4. The variance of the drag conductivity fluctuations (squares) plotted against 
temperature. The solid and dashed lines are calculated using Eq.[2]with the diffusive 
and ballistic asymptotes of r v , respectively Inset: r v extracted from the correlation 
magnetic field of the single-layer fluctuations, plotted against temperature. 
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